Draft version June 5, 2009 

Preprint typeset using I^T^X style cmulatcapj v. 2/19/04 



RADIUS DETERMINATION OF SOLAR-TYPE STARS USING ASTEROSEISMOLOGY: 
WHAT TO EXPECT FROM THE KEPLER MISSION 

Dennis Stello 1 ' 2 William J. Chaplin 3 , Hans Bruntt 1 , Orlagh L. Creevey 4 , Antonio Garci'a-Hernandez 5 , Mario 



o 
o 

(N 



Or 
I < 

o ■ 

c3 



> 

l> 
O 

o 

OS 

o 



X 



J. P. F. 



G. Monteiro 6 ' 7 , Andres Moya 5 , Pierre-Olivier Quirion 2 , Sergio G. Sousa 6,7 , Juan-Carlos Suarez 5 , 
Thierry Appourchaux 8 , Torben Arentoft 2 , Jerome Ballot 910 , Timothy R. Bedding 1 , 
Jorgen Christensen-Dalsgaard 2 , Yvonne Elsworth 3 , Stephen T. Fletcher 11 , Rafael A. Garcia 12 , 

GUNTER HOUDEK 13 ' 14 , SEBASTIAN J. JlMENEZ-REYES 4 , HANS KjELDSEN 2 , ROGER NEW 10 , CLARA ReGULO 4,15 , DAVID 



Salabert 



Thierry Toutain 



Draft version June 5, 2009 

ABSTRACT 

For distant stars, as observed by the NASA Kepler satellite, parallax information is currently of fairly low quality 
and is not complete. This limits the precision with which the absolute sizes of the stars and their potential transiting 
planets can be determined by traditional methods. Asteroseismology will be used to aid the radius determination of 
stars observed during NASA's Kepler mission. We report on the recent asteroFLAG hare-and-hounds Exercise#2, 
where a group of 'hares' simulated data of F-K main-sequence stars that a group of 'hounds' sought to analyze, aimed 
at determining the stellar radii. We investigated stars in the range 9 < V < 15, both with and without parallaxes. 
We further test different uncertainties in T e g, and compare results with and without using asteroseismic constraints. 
Based on the asteroseismic large frequency spacing, obtained from simulations of 4-year time series data from the 
Kepler mission, we demonstrate that the stellar radii can be correctly and precisely determined, when combined with 
traditional stellar parameters from the Kepler Input Catalogue. The radii found by the various methods used by each 
independent hound generally agree with the true values of the artificial stars to within 3%, when the large frequency 
spacing is used. This is 5-10 times better than the results where seismology is not applied. These results give strong 
confidence that radius estimation can be performed to better than 3% for solar-like stars using automatic pipeline 
reduction. Even when the stellar distance and luminosity are unknown we can obtain the same level of agreement. 
Given the uncertainties used for this exercise we find that the input logg and parallax do not help to constrain the 
radius, and that T c g and metallicity are the only parameters we need in addition to the large frequency spacing. It is 
the uncertainty in the metallicity that dominates the uncertainty in the radius. 

Subject headings: stars: fundamental parameters — stars: oscillations — stars: interiors 
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1. INTRODUCTION 

With the recent successful launches of the CoRoT 
(Baglin et al. 2006) and Kepler missions (Borucki et al. 
2008) we have entered a new era with strong synergy be- 
tween the fields of transiting exoplanets and stellar oscil- 
lations (asteroseismology). This synergy exists because 
of the common requirements for long uninterrupted high- 
precision time-series photometry, enabling the same data 
to be used for both purposes, providing complementary 
information. The investigation of stellar oscillations, and 
especially solar-like oscillations, provides a unique tool 
to probe the interiors of stars. In particular, asteroseis- 
mology can provide an independent radius estimate of a 
planet-hosting star, which can then be used to constrain 
the size of its transiting planet (s) (Christensen-Dalsgaard 
et al. 2007; Stello et al. 2007; Kjeldsen et al. 2009). 

The new space missions provide the first opportunity 
to measure solar-like oscillations in a large number of 
F, G and K main-sequence stars, which have photomet- 
ric amplitudes of only a few parts per million. These 
oscillations are global p modes, excited by near-surface 
convection. The frequency spectra (or p-mode spectra) 
of the oscillations show a characteristic spacing called 
the large frequency spacing, Ais, between modes of suc- 
cessive radial order n. This spacing scales as the square 
root of the mean stellar density and can therefore be used 
to constrain the stellar radius with very high precision. 
Investigating this potential for the Kepler mission has 
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been carried out in the framework of the asteroFLAG 
collaboration, whose aim is to develop and test robust 
tools for analyzing asteroseismic data on solar-like stars. 
AsteroFLAG includes members of the Kepler Asteroseis- 
mic Science Consortium (KASC; Christensen-Dalsgaard 
et al. 2007), and of the CoRoT asteroseismology team 
(Appourchaux et al. 2008). Much of the work within as- 
teroFLAG is founded on a series of hare-and-hound ex- 
ercises, where a group of 'hares' generates artificial data 
that the 'hounds' seek to analyze. 

The first part of the asteroFLAG investigation (Ex- 
ercise#l) reported by Chaplin et al. (2008) involved ex- 
tracting the large frequency spacing from artificial Kepler 
data. In this paper we build on those results to test how 
reliably we can obtain the stellar radius. We assume the 
availability of standard catalogue data (e.g. T e g, log g, 
V , and parallax), which will come from the Kepler Input 
Catalogue (Brown et al. 2005) or other sources. 

Of the more than 100,000 stars to be observed by the 
Kepler satellite, most will be observed at low cadence 
(30 min), but a few hundred targets will be observed in 
a high-cadence mode (1 min). The observing mode for a 
given target can be changed from low to high cadence 
if there are indications that the target shows transit- 
like events, which will benefit both the transit measure- 
ment as well as the supporting asteroseismic investiga- 
tion. While low cadence is sufficient to sample solar-like 
oscillations in evolved subgiants and red giant stars that 
have periods of hours to days, in this paper we focus 
on the prime targets of the Kepler mission - the main- 
sequence stars. These have oscillation periods of a few 
minutes to several tens of minutes, hence requiring data 
in the high cadence mode for a successful application of 
asteroseismology. 

§ 2 will provide a summary of the asteroFLAG hare- 
and- hounds Exercise#l, which feeds as input to Exer- 
cise^ reported in § 3. The discussion in § 3 includes the 
results. The detailed descriptions of the methods applied 
by each hound are given in § 4-8, followed by a general 
discussion in § 9 on how our results depend on the input 
parameters, and we highlight possible systematic errors. 
Finally, we give the conclusions in § 10. 

2. EXERCISE#1 RESULTS 

For convenience we give a short summary of the re- 
sults from the asteroFLAG Exercise#l (Chaplin et al. 
2008) relevant for this paper. That exercise concerned 
the extraction of the large frequency spacings from sim- 
ulated p-mode spectra of three artificial main-sequence 
stars designated Katrina, Boris, and Pancho. Their loca- 
tions in the H-R diagram are shown in Figure 1. For each 
star, 4-year time series were generated for all combina- 
tions of the following parameters: apparent magnitude 
V = 9, 11, 13, 15; inclination of rotation axis relative to 
the line of sight i = 0°,30°,60°; and rotation equal to 
one, two, and three times solar. In Figure 2 we show four 
examples of power spectra for Boris (V = 9, 11, 13, 15). 
While the oscillations and the characteristic frequency 
spacing, Av, are clearly seen for the brighter cases we can 
hardly sec the excess power for the fainter ones. How- 
ever, we saw from Chaplin et al. (2008) that the large 
frequency spacing can still be extracted in most cases in 
the tested magnitude range. The results from the exer- 
cise are shown in the last column of Table 1. We note 
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Fig. 1. — H-R diagram showing evolution tracks (gray lines) and 
the true position of the three artificial stars (star symbols), which 
were not known by the hounds. The la, 2a, and 3a error boxes 
show the positions of the stars according to the artificial "observed" 
traditional catalogue data, which were known by the hounds. For 
clarity we only show error boxes for cases K2, B2, and P2 (see 
Table 1). For details on how the catalogue data were generated see 
§ 3.1. 
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Fig. 2. — Fourier spectra of 4-year Kepler time series of the 
artificial star Boris. Each panel corresponds to different appar- 
ent magnitudes, and shows co-added spectra from subdividing the 
time series in four day subsets. The large frequency spacing, Au, 
between two modes of successive order n is indicated in the top 
panel. The catalogue data for each case are shown in Table 1. 
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that the ability to measure Av for the faintest case of 
Boris (V = 15), and perhaps V = 13, should be viewed 
with some caution. This is because the hounds in Ex- 
erciser had access to information about the frequency 
range in which they should measure Av from their anal- 
ysis of the brighter cases where power excess is clearly 
visible. Hence, in a true blind case, as with real data, it 
might be harder to establish Av for such faint cases of 
'Boris-like' stars. For the three faintest cases of Katrina 
(V = 11,13,15) the large frequency spacing could not 
be determined, preventing an asteroseismic constraint on 
the radius. 

The time-series data used in Exercise^l were gen- 
erated with the asteroFLAG simulator (Chaplin et al. 
2009, in preparation). Each time series included pho- 
tometric perturbations from p-modes, granulation, ac- 
tivity, photon noise, and instrumental noise. The stellar 
models were generated using the Aarhus stellar evolution 
code, ASTEC (Christensen-Dalsgaard 2008b) using the 
simple but fast EFF equation of state (Eggleton et al. 
1973) along with the Grevesse & Noels (1993) solar mix- 
ture using the OPAL tables at high temperature and the 
Kurucz (1991) table at lower temperature. The exact 
model values for L, T e g, and R are listed in Table 2. 
The frequencies of the p-mode signal were calculated us- 
ing the adiabatic pulsation code ADIPLS (Christensen- 
Dalsgaard 2008a). The amplitude and damping rates 
were from semi analytical models (Chaplin et al. 2005). 
We refer to Chaplin et al. (2008) for further details on 
the simulations and Exercise#l results. 

3. ASTEROSEISMIC RADIUS DETERMINATION: 
EXERCISE#2 

The aim of Exercise^2 is to use the large frequency 
spacings from Exercise^l, together with 'traditional' 
stellar parameters, to estimate the radii of the stars. 
While individual mode frequencies, and hence detailed 
asteroseismic analysis, are expected for the brighter 
stars, in this paper we only explore the benefit for the 
radius estimate from having access to the large frequency 
spacing. This work is an effort towards the development 
of robust algorithms for estimating radii of stars of differ- 
ent properties and noise levels in an automated pipeline 
fashion. It is the intention that the results are used to 
indicate what radius precision we can obtain on stars in 
the V — 9 — 15 range, as will be observed by Kepler. 

3.1. Stellar parameters 

In addition to the large frequency spacing obtained 
from asteroseismology by the first group of hounds (Ex- 
ercise#l), the second group of hounds (Exercise#2 - this 
paper) were given a set of traditional catalogue data for 
each star. These are given in Table 1. Each parameter, 
log 5, \og(Z/X), T c ff, and tt was obtained by adding to 
the true value (which the hounds did not know), a ran- 
dom number drawn from a normal distribution having 
zero mean and standard deviation equal to the adopted 
uncertainty on the parameter. Based on the traditional 
catalogue data we have estimated the locations of the 
stars in the H-R diagram, which are shown in Figure 1 
together with the true values (star symbols). In real- 
ity most of the atmospheric parameters, such as T e ff, 
log<7, and metallicity, will come from the Kepler In- 
put Catalogue, which is based on calibrated photometry. 



However, for many of the most interesting targets the 
photometric information will be supplemented by high- 
resolution spectroscopy (Latham et al. 2005). Hence, for 
each star at each magnitude we analyzed two cases. In 
one, the uncertainty in T e g was assumed to be 200 K, 
which is representative for the relevant stars from the Ke- 
pler Input Catalogue (T. M. Brown, priv. comm.). The 
other case assumed a smaller uncertainty in T e g , which 
would be the case if additional ground-based data, such 
as high- resolution spectra (Sousa et al. 2008), were avail- 
able. 

For the parallaxes we adopted Hipparcos-like preci- 
sions, with some extrapolation at the faint end to repre- 
sent a conservative estimate of the expected astrometry 
from the Kepler data. For V > 13 the Hipparcos-like 
parallaxes had large uncertainties, and we did not even 
generate predicted parallax data for V = 15. For this 
faint end, the luminosity could only be estimated very 
roughly using the traditional stellar parameters. 

The large frequency spacings quoted in Table 1 are the 
averages estimated from Exercise#l, and the uncertain- 
ties represent the scatter of the estimates. These values 
were made over results from all hounds on all datasets 
at each V value. We recall that datasets in Exercise^l 
were made for different internal rates of rotation, and 
different angles of inclination. Hence, the uncertainties 
given in Table 1 not only reflect the impact of reduction 
noise due to difference in analysis between hounds but 
also more subtle contributions from the different rota- 
tion and inclination. 

For the faintest case of Pancho (V = 15) there were 
outliers in the individual Av estimates returned by each 
hound in Exercise#l, and we decided to examine two 
versions of this star: one where Av was the mean of all 
the values returned by each hound, and the other where 
outliers were removed before calculating the mean. 

Finally, we also considered cases for all stars where 
the large frequency spacing was not used by the hounds. 
These are listed in the bottom part of Table 1, below the 
line. This allowed us to compare the radius estimates 
with and without the asteroseismic constraint. 

3.2. Analysis and results 

Each hound in Exercisc^=2 worked independently and 
used different methods to determine the stellar radii. De- 
tails are given in § 4-8. All methods rely on matching 
various parameters (e.g. L, Z, T c g-, and Av) from stellar 
models with the corresponding observables (e.g. V, log g, 
\og(Z/X), T c ff, 7r, and Av) given in Table 1 to estimate 
the radius. 

In Table 2 we list the radius estimates returned by the 
hounds as well as the model values. Not all hounds ex- 
amined all stars. By comparing the results from the top 
of the table (cases Kl-2, Bl-8, and Pl-10) with those 
at the bottom (cases K3-4, B9-16, and Pll-20) it is 
quite evident that including the large frequency spacing 
gives a significant improvement in our ability to deter- 
mine the radius. In most cases, including this asteroseis- 
mic constraint reduces the uncertainty in the radius to 
a few percent, which is an improvement by a factor of 
5-10. We note that in the bottom part of Table 2, the 
majority of the uncertainties from Hound-2 and a few 
from Hound-3, seem slightly underestimated compared 
to the deviations, (i? h ound - #modei)/-Rmodei, from the 
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TABLE 1 

Stellar input parameters given to each hound 
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true value. 

We summarize our main results in Figure 3. This 
shows the results listed in Table 2 for the stars that most 
hounds have examined, and which included the large fre- 
quency spacings. The uncertainties found by the hounds 
of roughly 1-3% generally seem to be consistent with the 
deviations that we see from the true values (see Fig. 3). 
We note a relatively small improvement in the radius es- 
timates when the uncertainty in T e s is significantly lower 
(left panels) than our base assumption of 200 K (right 
panels). This indicates that, provided we have the as- 
teroseismic constraint (the large spacing) , improving the 
estimate of T e ff is not particular important. In the fol- 
lowing sections we discuss in more detail the methods 
adopted by the six hounds, and refer to § 9 for further 
discussion on how our results depend on the input pa- 
rameters and their assumed uncertainties. 

4. RADIUS: RAPID ALGORITHM FOR DIAMETER 
IDENTIFICATION OF UNCLASSIFIED STARS 
(HOUND-1) 



In this section we describe the RADIUS pipeline devel- 
oped and utilized by DS to estimate stellar radii. The 
philosophy behind the development of this pipeline was 
to make it fast, robust and simple. It follows the same 
basic principle as the other methods described in the fol- 
lowing sections, namely comparing a number of observ- 
ables of a star with a grid of stellar models. All stel- 
lar models that are within 3cr of the observations in all 
four parameters T c g, L, Z, and Ais simultaneously are 
treated as being equally likely. Hence, a model either fits 
(is within 3a) and is accepted or it does not fit, in which 
case it is discarded. 

4.1. Models 

The selection of models (step (2) below) was based on a 
model grid that covered uniformly the entire range (±3er) 
in T e g, L and Z/X spanned by each star (see Table 1). 
The resolution in \og{Z/X) was 0.1 dex, corresponding to 
°~iog{z/x)i an d the resolution in mass was 0.01 Mq. The 
grid was generated with the Aarhus stellar evolution code 
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TABLE 2 

Results from the hounds 
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ASTEC using the simple but fast EFF equation of state 
(Egglcton ct al. 1973), a fixed mixing- length parameter, 
a = 1.8, and an initial hydrogen abundance of X = 0.7. 
We restricted the parameter space by fixing a and X. 
However, the effect of changing these parameters was ex- 
plored by other hounds (see § 5, 7, and 8), who showed 
the effects to be quite small. The opacities were calcu- 
lated using the solar mixture of Grevesse & Noels (1993) 
and the opacity tables of Rogers & Iglesias (1995) and 
Kurucz (1991) (T < 10 000 K). Rotation, overshooting 
and diffusion were not included. We used the adiabatic 
pulsation code ADIPLS to calculate for each stellar 
model. 

4.2. Pipeline approach and results 
The pipeline took the following approach: 

1. It determined the location in the H-R diagram 
by calculating L (Fig. 4; black cross and dotted 
3-cr error box). If the parallax was not known, 
as in the V = 15 cases, the luminosity was esti- 
mated from L/Lq = (i?/R Q ) 2 (T off /T offQ ) 4 , where 



R/R Q = (g/g )/(Av/Av Q ) 2 - 

2. It found the stellar models that matched within 
±3cr in L, T c ff, Z/X, and Av (see Fig. 4; colored 
symbols). 

3. Among all matching models, the two with the 
largest and the smallest radii, i? max and i? m i n , were 
identified. Lines of constant radius with these val- 
ues are shown in Figure 4 (dashed lines). 

4. The estimated radius was calculated as the average: 
R = (-Rmax, -Rmin), with uncertainty a R = (i? max - 
^min)/6 to accommodate that we found the two 
extreme models within ±3o\ 

In the following we give details for the two first steps: 

i To calculate L we interpolated the grid of bolomet- 
ric corrections by Lejeune et al. (1998), which takes 
the following input: T c g, log g, and [Fe/H] 17 . For 

17 Using bolometric corrections by Flower (1996), which only 
takes T c jj as input, did not change our final results significantly. 
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Fig. 3. — Radii of the artificial stars and their lcr-error bars estimated by each hound arc plotted in comparison with the true value 
indicated by the dashed lines. We show for each case the deviation in percent from the true value. Left panels show the results assuming 
<Tt ff = 200 K and right panels show more optimistic scenarios (or off (Katrina)= 100 K, cr Trff (Boris)= 25 K, and or off (Pancho)= 40 K.) 



that we first converted Z/X to [Fe/H] using the so- 
lar values Z Q = 0.0188 and X Q = 0.6937 from Cox 
(2000). We further adopted M bo i i0 = 4.746 from 
Lejeune et al. (1998). In the cases without parallax 
we used the solar values of T c g © = 5777 K, log g & = 
4.44 (Cox 2000) and A^ = 134.92 fiRz (Toutain 
& Frohlich 1992) to estimate L. 

ii To match models with observations we estimated 
Ai; by fitting to the radial modes of successive 



order n in the frequency interval 0.75-1.25 times 
i/ max (the frequency of highest power excess) , which 
we estimated from scaling the solar value (see § 9 
for a discussion on the best choice for the frequency 
range). The corresponding orders of n for each 
star were: 19-33 (Katrina), 15-27 (Boris), 13-22 
(Pancho). The calculation of Av was done for all 
models and the results stored in the grid, making 
the search for models that matched any given star 
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Fig. 4. — Hound- 1: Location of Boris (case B2) in the H-R di- 
agram. The stellar location (cross) and 3<r-error box (dotted) are 
shown in accord to the traditional catalogue data (see Table 1), 
while the true location of Boris is shown by the star symbol (not 
known by the hound). Colored symbols show the location of the 
stellar models with large frequency spacings that are equal to the 
observed value plus (squares) or minus (diamonds) 3<r. The dashed 
lines indicate lines of constant radius for the largest and smallest 
model that agree with observations within 3<r in all four param- 
eters: L, T e ff, Z, and Au. Note that models with fixed Z and 
Ai/ follow closely lines of constant radius. 



extremely fast. 

The results from this pipeline are summarised in Ta- 
ble 2 (column: Hound-1). 

5. SEEK: A STATISTICAL APPROACH FOR THE AARHUS 
KEPLER PIPELINE (HOUND-2) 

Here we describe the results obtained by a beta version 
of the seek routine developed by POQ for the Aarhus 
Kepler pipeline. Like radius (§ 4), seek uses an ex- 
tended grid of stellar models calculated with the Aarhus 
stellar evolution code (ASTEC) and the adiabatic pul- 
sation code ADIPLS to estimate the radius of stars. 

Like the RADIUS procedure, the seek routine was de- 
veloped for robustness, speed and simplicity, but we 
adopted a slightly more sophisticated approach. Instead 
of the brute force approach of RADIUS, which simply se- 
lects all models within 3er and calculates the radius as a 
simple average of the two extreme models (Fig. 4) , seek 
selects the set of best-fitting models based on a x 2 for- 
malism and estimates the radius by fitting a Gaussian to 
the distribution of radii of those models. 

5.1. Models 

The core of seek is the grid of models used to fit the 
observations. For this exercise we calculated and merged 
two regularly spaced subgrids, with various values for the 
mixing-length parameter and metallicity. The first sub- 
grid comprised 20 sets of evolutionary tracks, each with a 
different combination of the mixing-length parameter in 
the range 1.2 < a < 3.0 in steps of 0.6, and metallicity in 
the range 0.01 < Z < 0.03, in steps of 0.005. The second 
subgrid included 12 sets of tracks having 1.5 < a < 2.7, 



and 0.0125 < Z < 0.0275, both with the same resolution 
as in the first grid. Each of the 32 sets comprised 57 evo- 
lution tracks from 0.6 to 3.0 Mq. The spacing between 
the tracks was 0.02 M from 0.6 to 1.4 M Q and 0.1 M© 
from 1.4 to 3.0 M . 

In this beta version of the pipeline we used simple in- 
put physics to speed up the computation of the grid. In 
particular, the EFF equation of state (Eggleton et al., 
1973) was used. The final version will use the modern 
OPAL equation of state, covering a more extended re- 
gion of the metallicity and mixing-length parameter do- 
mains and include a range of hydrogen fractions, which 
in the current version was fixed at X = 0.7032. The mix- 
ture of elements was taken from Grevesse & Noels (1993) 
and opacity tables were from the OPAL package supple- 
mented by Ferguson et al. (2005) for low temperatures. 

5.2. Method and results 

We calculated the x 2 distribution of each star within 
our grid of models using: 

i— 1 ^ ' 

where Oi are the following five observed quantities: 7r, 
T e ff, logg, \og{Z/X), and Av, each with uncertainty <7j 
(see Table 1). The corresponding model quantities are 
Mj. In contrast to RADIUS, seek does not estimate the 
observed luminosity to compare with the models, but in- 
stead calculates a model parallax to match the observed 
apparent magnitude of the star, which is assumed to 
be exact. The color transformation is done using the 
VandenBerg & Clem (2003) tables. In the faint cases 
(V = 15) we had no parallax information, hence the par- 
allax was not included in the x 2 sum. We derived the 
large frequency spacing for each model by fitting to suc- 
cessive radial overtones of order n ^ 15 calculated with 
ADIPLS. These modes are in the asymptotic regime 
(Tassoul 1980), meaning they are close to being equally 
spaced, and hence well suited for the computation of the 
large spacing (see § 9 for a discussion on the best choice 
for the frequency range of the modes for the Ai/ calcula- 
tion). 

A typical result of the x 2 calculation is shown in Fig- 
ure 5 for Pancho (V — 9). This illustrates a slice of the 
multidimensional parameter space with a fixed metal- 
licity, Z — 0.0175, and a fixed mixing-length parame- 
ter, a — 1.5. It shows clearly a valley of best-fitting 
models along a contour of constant radius (in this case 
R ~ 1.7-Rq). While models with quite a range of masses 
and ages fit well to the observations, they all have roughly 
the same radius. However, due to a relatively flat bot- 
tom of the valley near the minimum there is no single 
best-fitting model. In other words the x 2 solution is de- 
generate. In addition, our grid has limited resolution and 
the ^-minimization problem is highly nonlinear. Hence 
we cannot assume that the model with the lowest x 2 m 
the grid is the absolute minimum of the problem. In- 
stead, we estimated the most likely value of the radius 
from a sample of best-fitting models for which x 2 was less 
than a given threshold, x 2 h . We did this by first split- 
ting the best-fitting models into bins of 0.005 Rq. For 
each bin we derived the total time spent by the models 
at that radius, which we denoted t(R). This gave us the 
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log y' = 1.175 1.979 3.763 3.587 4.392 5.196 6.000 



a = 1.5 Z = 0.0175 
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Fig. 5. — Hound-2: x' 2 calculated for Pancho (case P2) in the 
mass-age space for models with Z = 0.0175 and a = 1.5. Radius 
contours are shown as solid lines (solar units). The color follows 
the log x 2 as atop of the figure, r is the age normalized by the age 
at the giant branch, thus stars evolve parallel to the ordinate from 
the Zero Age Main Sequence at r = to the giant branch at r = 1. 
The kinks in the radius contours around 0.7 < r ^0.9 correspond 
to the blue hook in the H-R diagram. The dark red region in the 
upper right corner marks stars that are older than 14 Gyr. 

distribution, t(i?), of best-fitting radii, taking into ac- 
count how likely it was to have the observed star at each 
particular radius. 

If the threshold, Xth> ^ s well-chosen and the grid of 
models is dense enough to allow a large number of models 
to be within that threshold, t(R) will be a normal distri- 
bution centred around the most likely value of the radius. 
We obtained the best value for Xth by a n iterative pro- 
cess that minimized the difference between the median 
of the distribution, t(R), and the centre of a Gaussian 
fit. In other words, we minimized the fitting error on the 
free parameter i?o, called ARq, by fitting the following 
Gaussian to the distribution t(R) (see Fig. 6): 

P(R;R a ,*)= C M- iR ~f )2 h (2) 

where Rq is our solution for the radius. We adopted a 
as the uncertainty on this radius, which of course was 
much larger than the fitting error, Ai?o (see Press et al. 
1992, for details on the fitting procedure). We note that 
our result is not highly sensitive to the choice of xth as 
long as it is within a reasonable range of the best value. 
This is because the valley is relatively flat and has steep 
sides. However, if xih is chosen excessively big or small it 
will change the width of the distribution, t(R), and hence 
over- or underestimate the uncertainty on the radii. 

The seek routine is able to establish the radius of the 
stars without requiring an intensive computation for ev- 
ery new star to be analyzed (a few seconds per star on 
an Intel Pentium D machine). Our results are presented 
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Fig. 6. — Hound-2: Representation of the automatic SEEK 
process used to find stellar radii (Pancho; case P2). Left panel: 
Gaussian fit (red dashed line) to the distribution histogram t(R) 
(black solid line), of the radii of the best-fitting models, which 
fulfil x 2 < Xth = ^06. The distribution has been normalized to 
unity. Right panel: Fitting error on Ro as a function of the chosen 
threshold for selecting the best-fitting models. The red dot shows 
the threshold, x?y, = 306 that results in the best fit. 



in the Hound-2 column of Table 2. 

5.3. Future development 

If a given observable is not known, the extent of the 
grid will artificially set the minimum and maximum value 
of this observable. To take that into account, we will 
build our final grid to cover all reasonable values in metal- 
licity, mixing length, and mass expected for solar-like 
stars. This constrained parameter space can be used to 
extract information about a star for which the only well 
known parameter is Av. If we take a star with an un- 
certainty of 1 /iHz or less on Av and only assume that 
the star is located somewhere in the parameter space, we 
get a typical uncertainty of 5-10 % on the radius. This 
uncertainty is an absolute maximum if the large spacing 
is known. This capacity of providing a fast and reliable 
answer with minimum knowledge of the star is a powerful 
feature of the seek procedure. 

6. THE SHOTGUN METHOD (HOUND-3) 

In this section we describe the automatic pipeline de- 
veloped by HB, called SHOTGUN, to perform radius esti- 
mation of stars with or without asteroseismic constraints. 
Similar to radius and seek, it uses a grid of stellar mod- 
els, but differs by selecting a random sample of those 
models from which the radius is determined. While 
SHOTGUN is still quite fast (~ 1 minute per star), it is 
somewhat slower than RADIUS and seek. We note that 
the SHOTGUN method is the only one presented here that 
does not make use of stellar pulsation calculations, but 
obtains the large spacing from scaling the solar value. 

6.1. Models 
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We applied canonical scaled-solar BaSTI isochrones 
(Pictrinfcrni ct al. 2004) in version 4.1.0 with a mass- 
loss parameter r\ = 0.2. The grid we used included 
stellar masses from 0.5 to 5M Q with steps of typically 
0.02 M Q and metallicities 19 from Z = 0.001 to 0.04. 
We converted the input parameter \og(Z/X) to val- 
ues of Z by applying the hydrogen content of the Sun, 
X Q = 0.7395, in agreement with Y & = 0.2485 (Basu 
& Antia 2004) and Z Q = 0.012 (Grevesse et al. 2007). 
Each model had a set of values (age, Z, M, L, T c s) and for 
use in the following we calculated in addition R/Rq — 
(Teff/Teff^) -2 (L/L ) - 5 and the large spacing by scal- 
ing from the Sun, Ais = Av Q {M/M & )°- 5 (R/Rq)- 1 - 5 
(Kjeldsen & Bedding 1995). For Av Q we used 134.8 ^iHz 
from Kjeldsen ct al. (2008). 

6.2. Method and results 

To select an unbiased sample of models from the dis- 
cretely distributed BaSTI grid, shotgun first generates 
points from a non-discrete random distribution in the 
three-dimensional space L, T e g, and A^, and finds the 
corresponding best matching models. The software is 
an IDL code that generates a 'shot' that is comprised 
of 50 randomly chosen three-dimensional points or 'pel- 
lets', Pi — (L i} T c ff : i, Az/j), i G [1;50], with mean values 
and l-a Gaussian random errors as listed in Table 1. The 
value of Li is calculated from the V magnitude and the 
parallax, using bolometric corrections from Bessell et al. 
(1998), which depend mostly on T e g and only weakly 
on log g. When no parallax is available Li is ignored in 
equation (3) below. For each pellet the closest matching 
model in the BaSTI isochrones is found by identifying 
the highest value of the following weight function: 

Wi = 

OVT1 I" (Li-L B ) 2 (T cff ,,-T cffiB ) 2 (Ai/j-Ai/fl) 2 "Ln 
eX P 2al 2^ 2^ > 6 > 

where B are the models in BaSTI. This is similar to 
finding the lowest x 2 . We illustrate this in Figure 7, 
which shows an H-R diagram with BaSTI isochrones for 
metallicity Z = 0.03. The large square and 3-cr error 
ellipse represent Boris (case B2). The diagonal black 
line connects models on the BaSTI isochrones with the 
input large spacing of 135.9 /iHz. The solid circles are 
the 50 pellets generated by SHOTGUN and the open cir- 
cles mark were we find the highest value of the weights, 
Wi, for each pellet, hence representing the correspond- 
ing matched models. To illustrate all three parameters 
in Eq. 3 we show a similar diagram with the ordinate ex- 
changed with the large spacing. It is clear that Aiy dom- 
inates the weight function in this example because the 
luminosity is very uncertain, which explains why all 50 
pellets get matched to only 16 models within a narrow 
range apparently independent of the luminosity of the 
fired pellets. The output radius from SHOTGUN is calcu- 
lated as the mean value for the open circles in the BaSTI 
model grid while taking into account multiple hits on the 
same model, and the uncertainty is the l-cr RMS scatter. 
One could in principle include a term in Eq. 3 for the 

18 Available from http : //albione . oa-teramo . inaf . it/ 

19 Grid values were Z = 0.001, 0.002, 0.004, 0.008, 0.01, 0.0198, 
0.03, and 0.04. 









~- r 


V 


ibip 


r 134 




'-_ 3i 136 






- 138 






: 3. 769 3.762 3.756 
■ Z=0.020 log T«< 
c , , , 1 , , , 1 






3.80 3.78 3.76 3.74 

log T.„ 

Fig. 7. — Hound-3: H-R diagram showing isochrones from BaSTI 
for Z = 0.03. The square symbol marks the input value for Boris 
(case B2) and the black dots are the 50 pellets generated by the 
SHOTGUN software. The open circles mark the best model match 
for each pellet, denned as the highest value of the weights from 
Eq. 3. The 3-cr error ellipse is indicated. The inset shows a similar 
diagram with Au vs. T e g. 



metallicity, but the BaSTI grid is too coarse in Z. In- 
stead, we did the analysis for the two values of Z in the 
grid that bracket the input value. Interpolating between 
the two values gave the final result for the radius, includ- 
ing the uncertainty associated with the uncertainty of Z. 
Our results are listed in Table 2 (Hound-3). 

7. A GLOBAL FITTING APPROACH (HOUND-4) 

We now proceed to describe the method used by OLC 
to estimate the stellar radii of Katrina, Boris, and Pan- 
cho. Like the seek procedure (§ 5), this approach is 
formulated as a x 2 -minimization problem where the ob- 
servations are used to find the best stellar model that 
minimizes equation (1). However, this approach differs 
by only using four observations in the ^-minimization 
(A^, T ff, \og(Z/X), and logg), while the apparent mag- 
nitude, V, and the parallax measurement, tt, were used 
exclusively to estimate the mass range of the star. 

7.1. Models 

As is the case for seek and RADIUS, the stellar mod- 
els were derived using ASTEC and ADIPLS. The 
physics of the models included the OPAL95 opacity ta- 
bles (Rogers et al. 1996) the mixing-length treatment of 
convection (Bohm-Vitense 1958), with a fixed convec- 
tive core overshoot parameter dov = 0.25, but variable 
mixing-length parameter, a, nuclear energy generation 
rates, and the EFF equation of state (Eggleton et al. 
1973). 

7.2. Method 

Unlike the previously described methods, we did not 
use a pre-calculated grid of models. Instead, the % 2 min- 
imization was performed on one model at a time where 
derivative information was used to calculate the subse- 
quent models that minimized the x 2 , and so the only 
limitations in resolution were set by the output files of 
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the models (e.g. five decimal places in solar radius). It 
has the advantage of not being limited by the parameter 
space of a grid including the flexibility of changing the 
input physics of the models. The disadvantage is that 
it needs to be run for various initial guesses to ensure a 
global minimum is found, with each run requiring new 
models to be calculated, which is more time consuming 
than the grid approach. 

Because we only had four observations to estimate the 
five stellar model parameters (mass, age, A, Z, and a), 
we had to fix one of the parameters, and we chose this 
to be X. In order to also explore the parameter X, the 
minimization was repeated while the fixed X took values 
between 0.68 and 0.74 in steps of 0.01. 

We use the Levenberg-Marquardt algorithm, along 
with singular value decomposition for the minimization 
(Press et al. 1992). The minimization is initiated giv- 
ing each observed parameter, Oi, their uncertainties, Oj, 
and an initial guess of the model parameters (mass, age, 
X, Z, and a) as input. The solution is found when x 2 
has reached 0.5 or after four iterations, which is suffi- 
cient for a stable solution. The frequency range used to 
derive Ai> of the model depended on the star and the 
magnitude. Examples for the V = 9 cases are Katrina 
6058-7915 //Hz, Boris 2500-4100 //Hz, and Pancho 1400- 
1978 //Hz. 

To avoid local minima problems, various runs were exe- 
cuted using different initial guesses of the mass, the range 
of which we estimate by deriving My from V and w and 
subsequently using the relations between My and mass 
from Allen (1973). The scheme automatically found the 
stellar model that minimized x 2 f° r each initial guess. 
Because we had various runs of the minimization with 
different initial masses (Nm) and different initial X val- 
ues (Nx), we essentially ran Nm x Nx minimizations 
and thus obtained this same number of best-fitting mod- 
els (see Figs 8 and 9). We determined the radius of each 
of the best-fitting models and chose the radius values 
from the models whose x 2 is below a suitable threshold. 
These models correspond to the filled circles in Figures 8 
and 9. The quoted radius in Table 2 is the mean value 
of these, while the uncertainty is defined as half of the 
difference between the largest and smallest radius value. 

7.3. Results 

Figure 8 shows examples of the masses and ages of the 
best-fitting models for various runs for Pancho (V=9). 
The initial guesses of the mass were 1.25, 1.30, 1.35, and 
1.40 M©, while the X value was also changed for each 
run. There are a total of 4 x 7 minimizations, and there- 
fore we have a total of 28 best-fitting models, each rep- 
resented by a circle whose size is inversely proportional 
to the x 2 value as defined by equation (1) (i.e. a larger 
point means a better fit). The filled circles are those 
models that give a x 2 < 3.9 2 . For each (almost verti- 
cal) ridge in Figure 8, X varies between 0.68 and 0.74. 
The dotted lines connects the results for X = 0.69, 0.71 
and 0.73. By inspecting the figure, it is clear that there 
are correlations among the parameters, meaning that we 
can obtain the same x 2 value (same dot size in Fig. 8) 
for a range of parameter combinations. The filled circles 
are those models that we accept as the best models, but 
still mass, age and X span large parameter ranges. This 
clearly demonstrates that some of the global parameters 
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Fig. 8. — Hound-4: Best-fitting parameters of mass and age for 
various minimization runs for Pancho (case P2). The larger circles 
represent the models with the lower \ 2 values. The filled circles 
represent the models whose \ 2 < 3.9 2 . The dotted lines connect 
the results for selected values of X. 
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Fig. 9. — Hound-4: Radius and effective temperature for the 
same models as in Figure 8. The dotted line is the observed T e g. 



cannot be significantly constrained by the observations 
that we have available in this exercise. For example, we 
can see the age of the best-fitting models varies between 
a significant range of 1.4 and 2.1 Gyr, and the mass varies 
from about 1.30 to 1.45 M Q . In fact, the 1-a uncertain- 
ties on the model with the lowest x 2 value are er(Mass) 
= 5% , a(age) - 21%, a(Z) = 23% and a(a) = 10%. 

In Figure 9 we show the corresponding radii and effec- 
tive temperatures of the same models as shown in Fig- 
ure 8. The dotted line indicates the observed T e g. All 
of the models fall within ±3<7T eff • In agreement with the 
results we found using seek and RADIUS most of these 
models have roughly the same radius despite spanning a 
large parameter range in mass, X, and age. 

This automatic minimization scheme was repeated for 
each of the stars and for each V value, and the results 
are listed in Table 2 (column Hound-4). 

8. NON AUTOMATED APPROACHES 

While the four methods described in the previous sec- 
tions (§ 4-7) were all automated, three relying on large 
existing grids of stellar models, we now turn to meth- 
ods that involve generating small sets of models covering 
more limited ranges of the parameter space, and itera- 
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tively narrowing down the parameter space of best-fitting 
models manually. The fundamental idea - fitting model 
parameters to observations - is still the same, but instead 
of being developed to the level of fast pipeline algorithms, 
which can be applied routinely on large sets of stars, the 
following methods again serve to demonstrate the abil- 
ity to determine precise radii. In addition, these meth- 
ods use different stellar evolution and pulsation codes to 
those in § 4-7, allowing a valuable comparison. 

8.1. The Granada approach (Hound-5) 

In this section we present the methodology adopted by 
JCS, AM, and AGH to estimate the radii of the three 
artificial stars (see Table 1). In contrast to the other 
hounds, we did not use the apparent magnitude and the 
parallax information. We therefore have only one case for 
each star, for which we adopted the average large spacing 
and the largest uncertainty listed in Table 1 (case K2, B6, 
and P8). 

8.1.1. Modelling 

To model the stars we constructed standard non- 
rotating models from late pre-main-scqucncc to terminal- 
age main sequence (where core hydrogen fusion ceases). 
We note that this does not cover all evolution states that 
are consistent with the observed parameters and their 
error bars, especially for Pancho, and our results may 
therefore have systematic errors. The evolutionary stel- 
lar models were computed with the CESAM code (Morel 
1997). Opacity tables were taken from the OPAL pack- 
age (Iglesias & Rogers 1996), complemented at low tem- 
peratures (T < 10 4 K) by the tables by Alexander & 
Ferguson (1994). The atmosphere was constructed from 
a grey Eddington T — r relation (Eddington 1926). Con- 
vection was treated with a local mixing-length model 
(Bohm-Vitense 1958) and we assumed various values for 
the mixing-length parameter a = l/H p , where I is the 
mixing length and H p is the pressure scale height. In par- 
ticular, values between 0.5 and 2 were explored. For the 
overshoot parameter e?ov = 'ov /H p (low being the pen- 
etration length of the convective elements) we assumed 
values between 0.0 and 0.3. Such ranges of parameters 
cover the range of parameters generally adopted for solar- 
like stars. 

For the computations of adiabatic cigcnfrequcncies two 
oscillation codes were used: FILOU (Tran Minh & Leon 
1995; Suarez 2002; Suarez & Goupil 2008) and GraCo 
(Moya et al. 2004; Moya & Garrido 2008). 

8.1.2. Method 

First, we calculated about 30 initial models that cov- 
ered the observed ±3cr range of each parameter, logg, 
\og(Z/X), and T e s, while making sure they also covered 
various values of mixing length and overshoot. From 
those initial models we then selected 3-5 representa- 
tive models with various masses and evolution states for 
which the large spacing was also within 3er of the ob- 
served value. The exact number of representative models 
varied from star to star. We consider the large spacing 
as the strongest bound on the radius. Hence, in order to 
get a better radius estimate we made an additional 3- 
4 models around each representative model, which were 
then all within 3a of the observed log g, \og(Z/X), T e g, 



and Au, with the majority within la of the latter. This 
provided a total of about 9-20 representative models. 
Our final estimate of the radius (Table 2, column Hound- 
5) is a simple average of the radii of those represen- 
tative models (i? = (J? mo dei)) and the corresponding 
uncertainty is the maximum deviation from the mean 

(a R = MAX(i? m odcl - (iJmodel)))- 

8.2. The CAUP approach (Hound-6) 

In this section we describe the procedure used by SGS 
and MJPFGM to estimate stellar radii. This approach 
is similar to the one described in § 8.1, in that an initial 
set of stellar models was made based on the observed pa- 
rameters and their uncertainties. Then, by testing which 
models fit best (in the logg, L, T c g parameter space), 
further refinement of the range of parameter values were 
made by calculating additional models. Finally pulsation 
models were derived to further constrain the parameter 
space using the strong bound from Av. 

Of the 13 stellar cases listed in Table 1, the procedure 
described here was only applied on one star (Boris, case 
B2, B4, and B6). The method does not rely on a pipeline 
algorithm, but can be improved and used for all the stars. 
Unlike the methods described in the previous sections, 
this approach only considered stellar models to be valid 
if they were within one sigma of the observed logg, L, 
and T e ff simultaneously. 

8.2.1. Models 

To obtain the evolution tracks we used the CESAM 
code version 2K 20 (Morel 1997; Morel & Lebreton 2008). 
The Livermore radiative opacities were used (Iglesias & 
Rogers 1996) complemented at low temperatures (T < 
10 4 K) by atomic and molecular opacity tables from Ku- 
rucz (1991). The opacities were calculated with the so- 
lar mixture of Grevesse & Noels (1993), convection was 
described by the mixing length model (Bohm-Vitense 
1958), and overshooting (dov) was not used since this 
parameter is poorly known, in particular for stars below 
2 M© (Ribas et al. 2000). The nuclear reaction rates were 
from the Nuclear Astrophysics Compilation of REaction 
Rates (NACRE) (Angulo et al. 1999). For the equation 
of state the OPAL 2005 tables 21 were used. The atmo- 
sphere of the stellar models was based on a Hopft law 
(Mihalas 1978) where convection was not included, and 
radiative transfer was considered to be independent of 
the radiation frequency (i.e. the grey case was assumed). 
We derived frequencies of the models with the pulsation 
code POSC (Montciro 2008) 

8.2.2. Method 

By using the parameters in Table 1 we located the 
star and its 1-a error box in both the H-R and logg 
vs. T e ff diagrams. To locate the position in the H-R 
diagram we estimated the stellar luminosity, computed 
from the parallax and the apparent magnitude, using the 
calibration by Flower (1996) to obtain the bolometric 
correction. In the calculation we adopted the solar ab- 
solute magnitude of My = 4.81 (Bessell et al. 1998) and 
a bolometric correction for the Sun of 0.08 mag (Flower 

20 CESAM2K available at www.oca.eu/cesam/ 

21 Tables available at http : //phys . llnl . gov/Research/0PAL/E0S_2005/ 
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Fig. 10. — Hound-6: Radius vs. large spacing for the 188 best- 
fitting models within one sigma in log g, log L, and T c g (Boris, case 
B2). The vertical line shows the measured Av and the horizontal 
line is the corresponding radius found by this approach. 



1996). The uncertainty on the luminosity comes mainly 
from the parallax. We then built a grid of stellar mod- 
els, for different masses, M, initial helium content, Yq, 
and mixing length, a, while having the initial metallicity 
ratio Z /X fixed and equal to the observed Z/X. In 
the following we use Boris (case B2) as an example to 
describe the process of finding valid fitting models. 

In an iterative process, we calculated three grids of 
models (see Table 3; grid 1, 2, and 3) to determine 
the range in M, Y , and a, that made the models fall 
within ±1(T of the observed log 3, T c s, and L, simulta- 
neously. In this process we made evolution tracks with 
all the possible combinations of M, Yq, and a within 
each grid. We finally computed two more grids (4 and 5) 
with higher resolution but covering roughly the same pa- 
rameter space as the third grid. With our adopted fixed 
ratio \og(Zo/Xo) = —1.6, the range in Yq corresponded 
to Z = 0.017-0.019, and X = 0.68-0.75. We note that 
this is a relatively small change in metallicity compared 
to the uncertainty quoted in Table 1. Hence, it is possible 
that we do not see the full effect on our radius estimate 
from the uncertainty in the metallicity. 

From all the grids, roughly 300 models from 32 tracks 
were within ±lcr in log g, logL, and T e g simultaneously. 
From this selection of models we picked one reference 
model on each track, and compute the pulsation frequen- 
cies in order to obtain the large frequency spacings, Av, 
which wc derived by fitting successive radial overtones in 
the frequency range 2400 4000 //Hz. Depending on the 
models the corresponding radial orders fell in the range 
n = 16 30. To save time computing pulsations models, 
we estimated Av for the rest of the models along each 
track using the following scaling relation: 

Av = (M/M ref ) 1 / 2 (i?/ J R ref )-3/2 Az/ref) 

where Av rc {, M re f, i? rc f are the large frequency spac- 
ing, the mass, and the radius of the respective reference 
model. 

We finally selected the best two thirds of the models 
(188 models in the case B2), which showed the smallest 
deviation between our estimated large spacing and the 
observed value. 

In Figure 10 we plot the relation between the large 
spacing and the radius for those models. We made a 
simple linear fit to the points presented in Figure 10 and 
used it to estimate the radius of the star, by evaluating 
the value of the fit at the observed large spacing (see 
vertical and horizontal lines in Figure 10). 



To estimate the uncertainty on the radius we combined 
the uncertainty from the fit and the uncertainty in the 
observed large frequency spacing. We obtained the latter 
by multiplying the uncertainty on the observed large fre- 
quency spacing with the slope of the fit. The uncertainty 
from the other observables are included in the uncer- 
tainty of the fit, since the points used for the fitting were 
obtained considering all models within the error boxes. 
The results for the cases B2, B4, and B6 can be seen in 
Table 2 (column Hound-6). 

9. DISCUSSION 

In this section we will discuss how our radius estimates 
are constrained by our input parameters and their as- 
sociated uncertainties. We further discuss the possible 
systematic errors in the results. 

9.1. Dependencies on input parameters 

Figure 4, which is the output from the radius pipeline, 
illustrates to some extent the effects of the various input 
parameters. All models that agree with the observations 
within 3cr (colored symbols) arc bound by T e ff, Av, and 
the metallicity. The gravity and parallax (when avail- 
able) are in general so uncertain that they do not con- 
strain the luminosity enough to exclude any models not 
already excluded by T e s, Av, and the metallicity. Only 
in a few of the brightest cases (V — 9) does the paral- 
lax exclude some additional models, reducing an slightly. 
Hence, despite the increasing uncertainty in the luminos- 
ity for the fainter stars we basically get the same radius 
and an (see Table 2). It is in fact mostly the increase 
in (tai/ that is responsible for the slight increase in an 
towards fainter stars. The redundancy of log g is further 
confirmed by re-running the analysis for a few of the 
stars with increased uncertainty a\ ogg = 0.2, and again 
wc obtain the same results within the uncertainty. 

Again, if we use Figure 4 as illustration we see that <ta^ 
and <7T eff contribute roughly equally to the uncertainty 
in the radius. However, it is evident that the largest 
contribution comes from the uncertainty in the metal- 
licity. It is a log ( Z / X ) that is responsible for increasing 
the acceptable range of radius values. It will therefore 
require a much lower crT cff and av to reduce an by the 
same amount as we will achieve by a slightly better deter- 
mined metallicity. We regard the applied a\ g(z/x) =0.1 
as realistic, but if one can reduce the uncertainty in the 
observed metallicity by a factor of two to four, we can 
reduce the uncertainty in the radius down to <~ 1%, de- 
pending on the position of the star in the H-R diagram. 
Only in the case of a very poorly determined Av does 
at\ v contribute more than <7i og (z/x) to the final an (see 
Table 2; cases P7 and P8). 

We further note from Figure 4 that the models of con- 
stant Av for a fixed metallicity, say the black squares, 
run almost parallel to lines of constant radius (dashed) . 
This explains why decreasing aT e!! has so little effect on 
an (see Table 2; cases Bl, B3, B5, and B7 versus B2, 
B4, B6, and B8). This parallel alignment is even more 
pronounced for the hotter star Pancho, which therefore 
shows no decrease in an as we decrease 0T eff . For the 
cooler star Katrina, the alignment is not as good, and a 
better-determined T e g does improve the radius estimate. 
A similar illustration is shown in Figure 5, which shows 
that the best fitting models (low x 2 ) lie along contours 
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TABLE 3 
Stellar model grids (Hound-6) 



Grid 


1 2 3 4 5 


M/Mq 

Y 

a 


0.94 0.97 1.03 1.06 
0.24 0.26 0.28 0.30 
1.20 1.40 1.60 1.80 


1.02 1.03 1.04 
0.25 0.26 0.27 
1.50 1.60 1.70 


1.08 1.10 1.12 
0.24 0.26 0.28 0.30 
1.20 1.40 1.60 1.80 


1.11 1.12 1.13 1.14 
0.23 0.24 0.25 0.26 0.27 
1.30 1.40 1.50 1.60 


1.08 1.09 1.10 

0.23 0.24 0.25 0.26 

1.30 1.40 1.50 1.60 1.70 1.80 



of constant radius despite having a large range in age 
and mass. The tight correlation between R and Av seen 
in Figure 10 further illustrates that this relation is not 
strongly dependent on the initial physics (mass, mixing 
length, initial helium abundance) of the star, which un- 
derpins that Av is a powerful measure to constrain the 
stellar radius. 

In summary, provided we have access to the large spac- 
ing, only metallicity and to less extent T e g provide ad- 
ditional constraints to the radius, while the parallax and 
log g are redundant due to their large uncertainties. 

9.2. Systematic errors 

For a given star, the large spacing is not strictly the 
same at all frequencies in the p-mode spectrum, but 
varies smoothly with frequency. Hence, to correctly com- 
pare observations and models, it is necessary to use the 
same frequency range in both cases when Av is derived. 
However, this was not attainable in this study because 
the adopted observed large spacing was the average of 
all values returned by each hound of the previous Ex- 
crcisc^l, each of whom used different frequency ranges. 
Fortunately, it is straightforward to employ a common 
frequency range for the observations and the models in 
the final pipeline for analyzing real Kepler data. 

We further note that all hounds used only the ra- 
dial modes to derive Av from the models, while the 
observed values included information on all detectable 
modes. This introduced a systematic error on the radius 
measurement because modes of different angular degree 
I have slightly different spacings. For a solar model this 
difference is 0.5% at most. Since Av oc p 5 we can get 
a rough estimate of the corresponding systematic error 
on the radius of < 0.3%. However, this again can be 
corrected for and the effect be eliminated. 

When fitting stellar evolution tracks to the stellar loca- 
tion in the H-R diagram, a common strategy is to restrict 
the search to the l-a error box (L, T e ff). Our investiga- 
tion clearly shows that this approach can lead to wrong 
results. With the very precise measurements of Av, we 
saw stars where traditional parameter values (e.g. L, 
T e f[) would deviate by more than lcr, in order to match 
Av within 3<r, which explains why some of the stars fall 
outside their l-a error box in Figure 1 . Statistically this 
is of course not surprising, but we still find it important 
to stress the point. 

In addition to the statistical aspect, inconsistencies 
can arise if there are systematic differences between the 
physics included in the modeling of the hares and the 
hounds or between the real star and the models. The 
stellar model codes used by the hares and hounds have 
been developed over many years and are well tested. For 
example, these codes have been shown to give almost 
the same radii and luminosities, with differences on the 



order of 0.05% (Lebreton et al. 2008). Further, com- 
paring results from the various hounds suggest that the 
detailed input physics has a relatively small effect on the 
radius estimate. Recently, Kjeldsen et al. (2009) sug- 
gested that the well-known offset of roughly 1% between 
the observed Av of the Sun and the value derived from 
stellar models, which comes from improper modeling of 
the near-surface layers of the Sun, can be corrected for 
empirically. This offset, which is likely to affect mod- 
els of other stars as well, has not been accounted for in 
this investigation because it requires knowledge of indi- 
vidual mode frequencies. The systematic effect on the 
asteroseismic radius estimate from this offset is therefore 
expected to be roughly 0.7% for a Sun-like star. 

Finally, we note that the hounds were given log(Z/X) 
and their results therefore did not depend on the adopted 
solar value to any large degree. However, in reality one 
needs to convert the observed metallicity to X, Y, and Z 
to compare with the models, and this conversion depends 
on the assumed solar metallicity. The range of Zq values 
currently accepted, 0.012 < Zq < 0.017 (Grevesse et al. 
2007; Caffau et al. 2008) will correspond to the same 
range in Z for a star similar to the Sun. This change, 
AZ — 0.005, will change the estimated radius by about 
1-2%. 

10. CONCLUSIONS 

In an extensive hare-and-hounds exercise, we have an- 
alyzed artificial solar-type main-sequence stars, which 
are representative for targets observed during the Kepler 
mission. The stellar radii found by each hound, using 
different approaches and codes for stellar evolution and 
pulsation calculations, are in good agreement with one 
another and with the true values. 

We are able to determine the radius with a precision 
of a few percent, and are confident that radii can be ob- 
tained with < 3% accuracy for solar-like Kepler targets 
on a routine basis using automatic pipeline reduction. 
We demonstrate that the reason for this is the strong re- 
lation between the radius and the large frequency spac- 
ing, which can be measured to very high accuracy. Our 
results further show that the radius is only weakly depen- 
dent on the input physics of the stellar models, including 
the equation of state and the mixing-length parameter. 

The uncertainty in the radius determination is mostly 
dominated by the uncertainty in the stellar metallicity, 
which translates to an uncertainty in the stellar mass. 
In most cases the parallax and the surface gravity are 
too uncertain to constrain the radius, so in essence only 
the uncertainty in effective temperature and the large 
frequency spacing adds further to the final uncertainty 
in the radius. However, we conclude that we do not need 
to know T e ff to a very high precision in order to obtain a 
good estimate of the radius, provided we have access to 
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Because the Kepler mission will provide high accuracy 
data for many binary systems, there will be cases where 
we will obtain both an 'asteroseismic radius', and a 'pho- 
tometric radius' giving opportunities to independently 
test the accuracy of our radius estimates. 

While these results rely on 4-yr time series, similar 
results arc expected from time series of only 1-3 months 
but with a limiting magnitude of about V ~ 11. Details 
on limiting magnitudes for determination of the various 
frequency spacings is the aim of ongoing asteroFLAG 
exercises. 
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